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DATA  DESCRIPTION  AND  QUALITY  ASSESSMENT  OF 
IONOSPHERIC  ELECTRON  DENSITY  PROFILES  FOR  ARPA  MODELING  PROJECT 
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National  Geophysical  and  Solar-Terrestrial  Data  Center 
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ABSTRACT 

This  report  presents  a  description  of  the  automated  method  used  to  produce  electron  density 
(N(h))  profiles  from  ionograms  recorded  on  35mm  film  and  an  assessment  of  the  resulting  data 
base.   A  large  data  base  of  about  30,000  profiles  was  required  for  an  ionospheric  modeling 
project.   This  motivated  a  search  for  an  automated  method  of  producing  profiles.   The  auto- 
mated method  used  is  fully  described,  the  resulting  data  are  given  a  quality  grade,  and  the 
noon  and  midnight  profiles  are  presented.   Selected  portions  of  this  data  base  are  compared 
with  profiles  produced  by  the  standard  profiling  method  in  use  by  the  Environmental  Data 
Service  at  Boulder,  Colorado. 

I.   INTRODUCTION 

The  derivation  of  electron  density  distribution  profiles  from  film  ionograms  by  the  Paul  [1967] 
method  as  implemented  by  the  computer  program  of  Howe  and  McKinnis  [1967]  has  been  an  ongoing  program 
in  the  Department  of  Commerce  since  early  in  1960.  However,  scaling  the  ionograms  and  preparing  the 
data  have  in  the  past  been  a  very  expensive,  tedious  and  time-consuming  effort. 

In  1972,  there  was  a  requirement  for  a  large  quantity  of  profiles  for  an  ionospheric  modeling 
project.  This  motivated  a  search  for  a  more  efficient  data  processing  system  by  Lerfald  et  al.   [1974] 
and  his  associates  at  the  Space  Environment  Laboratory  of  NOAA.  The  technique  decided  on  was  (1)  pre- 
examination  and  expert  commentary  of  samples  of  the  ionograms  to  be  analyzed,  (2)  digitization  by  a 
nearly  automatic  laser  beam  digitizer  and  (3)  introduction  of  a  suitable  model  for  "unseen"  ionization. 
The  difficulties  inherent  in  this  and  any  other  automated  method  of  producing  electron  density  profiles 
are  how  to  get  the  meaningful  ionogram  information  accurately  into  digital  format,  and  what  assumptions 
are  appropriate  when  information  is  missing  on  the  ionogram.  Generally,  ionograms  reveal  little  or  no 
information  at  low  frequencies  (below  1.5  MHz),  and  in  "unseen"  regions  (valleys  and  topside)  or  in 
"seen"  regions  where  the  electron  density  gradient  is  nearly  zero.  Solutions  to  both  of  these  problems 
are  important  if  we  are  to  calculate  reasonably  accurate  "true"  heights. 

Ideally  the  inversion  process  would  start  at  zero  frequency  but  in  actual  practice  this  is  not 
possible.  In  the  evaluation  of  the  inversion  integral  for  every  probing  frequency,  there  is  a  singularity 
where  the  plasma  frequency  is  equal  to  the  probing  frequency  for  the  o  component.  This  singularity  has 
an  effect  which  depends  greatly  on  the  slope  (dN/dh)  at  that  point.  If  we  could  choose  the  correct 
height,  frequency  and  slope,  we  would  be  able  to  "start"  the  inversion  at  this  point  and  thereby  closely 
approximate  the  actual  profile  as  long  as  the  rest  of  the  profile  was  monotonic.  However,  there  are   few 
data  in  the  ionogram  that  can  give  guidance  to  these  choices. 

Many  different  models  (or  "starts")  can  and  have  been  used  in  an  attempt  to  account  for  the  "unseen" 
ionization.  For  the  purpose  of  the  joint  ARPA/N0AA  study,  a  diurnally  varying  D-,  E-,  and  lower  F-region 
model  was  used  and  a  monotonic  distribution  assumed  throughout  all  regions  of  the  ionosphere.  Thus  the 
major  part  of  this  publication  will  be  concerned  with  reporting  how  well  the  profiling  procedure  was 
automated  by  Lerfald  et  al.     An  analogous  hut  different  technique  was  adopted  by  Dr.  John  S.  Nisbet  and 
colleagues  at  Pennsylvania  State  University  for  mass  studies  for  a  different  time  period  in  connection 
with  the  same  modeling  project. 

The  ionospheric  modeling  project  for  which  these  large  quantities  of  N(h)  profiles  were  prepared 
was  suspended,  but  the  N(h)  data  have  been  deposited  in  World  Data  Center  A  for  Solar-Terrestrial  Physics. 
This  report  describes  these  data  and  summarizes  the  Lerfald  data  processing  and  analysis  system.  It  also 
describes  and  illustrates  a  data  quality  assessment  scheme  which  has  been  applied  to  the  archived  Lerfald 
data.  The  Appendix  gives  graphs  of  the  noon  and  midnight  profiles  from  the  Lerfald  data  set,  illustrative 
of  the  complete  and  detailed  data  available  from  World  Data  Center  A,  which  is  operated  by  and  is  col- 
located with  NOAA's  National  Geophysical  and  Solar-Terrestrial  Data  Center. 


2.  DATA  DESCRIPTION 

Ionograms   for  the  Regular  World  Days   (RWD)   in  March,  June,  September,  and  December  1964   (IQSY) 
were  originally  selected  for  producing  the  electron  density  distribution  data  base  necessary  for  the 
ARPA  ionospheric  modeling  project.     However,   due  to  the  high  level   of  geomagnetic  activity  during 
September  and  December,   the  period  of  coverage  was   extended  in  these  two  months   to  include  geomagnetical ly 
quiet  days.     Except  for  September  these  include  the  dates   recommended  by  the  XVII  General   Assembly  of  URSI 
(Aug.  1972)   for  data  reduction   in  Recommendation   URSI   Commission   III . 1-Profi les  of  Electron   Density. 
For  stations  with  ionosonde  sweep  schedules  of  one  each   15  min,   2,000  ionograms   per  station  would  have 
been  taken  during  these  periods.     The  distribution  and  magnitude  of  the  quiet  and  disturbed  times  as 
indicated  by  Ap  and  Kp  can  be  seen   in  Table  1.   A  survey  of  the  available  ionogram  data  during  these  peri- 
ods  revealed   15  stations    (shown  in  Table  2)  with  the   15-min  sweep  schedule  and  reliable  data.     These 
15  stations  operating   continuously  for  the  22  days  would  have  established  about  30,000  ionograms   available 
for  processing  in  this  modeling  project.     Only  18,226  interpolated  profiles  were  actually  produced, 
because   ionograms  were  lost  due   to  equipment  malfunctions  or  were  too  complex  to  interpret. 

The  profiles  are  archived  in  computer  compatible  form  at  the  World  Data  Center  A  for  Solar-Terres- 
trial  Physics  in  two  formats  on  magnetic  tape.     Examples  of  these  formats  are  shown  in  Tables  3  and  4. 
Table  3  gives  a  sample  printout  of  the  profile  data  in  terms  of  the  points  used  in  the  calculation. 
Table  4  is   a  sample  of  the  profile  interpolated  to  true  height  intervals  of  ten  kilometers  as   they 
appear  in  card  or  card  image  format.     Each  profile  is   contained  on  two  or  more  cards. 

Other  output  formats   can  be  provided   upon   request,   such  as  shown   in   Figure   1  and  Table  5. 

Table  1.        Solar- Geophysical  Data  for  Days  Scaled  in  1964 
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Mar 


Jun 


Sep 


Dec 


Day 

Daily 

Three-hr 

Gr. 

Intervals 

Kp 

No . 

Ap 

magnetic 

inde 

X 

1 

2 

3 
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6 

7 

8 
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77 

10 

2  + 

2o 

4- 

4- 

2o 

1+ 
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1- 
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78 

3 

0+ 

0  + 

2- 

1  + 

1- 

1- 
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0+ 
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79 

3 

1- 

0+ 

2- 

0+ 

lo 

1- 

1- 
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4 

2- 

lo 

1- 
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0+ 

lo 

lo 
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4 

0+ 
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2- 
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2- 

1- 

0+ 

0+ 
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1- 

0+ 
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3- 

18 
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2- 
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2- 

19 
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0+ 

1- 

1- 
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0+ 
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2U 
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2 

Oo 
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0o 

0+ 

0+ 
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1- 

1  + 

21 
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5 

0+ 

1- 

1+ 

2- 
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7- 

8- 

3+ 

2+ 
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2o 

2+ 
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2o 
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4- 
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9 
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Oo 

1- 

1- 

13 

348 

9 
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1+ 

1- 
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4- 

3+ 
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14 

349 
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2o 
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15* 

350 

8 

3o 

3+ 
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lo 

2- 

lo 

lo 

0+ 

16* 

351 

15 

3- 

2  + 

2o 

3- 

5+ 

2o 

2- 

3- 

17* 

352 

15 

5- 

4- 

3o 

3o 

1+ 

2- 

2o 

2+ 

''Regular  World  Days 


Table  2.      A  L 


zst  of  Ground  Stations  with  Geographic  and  Geomagnetic  Coordinates 


Station    Latitude 

Long 

i  tude 

Geomagnetic* 

Magnetic* 

Station         Code 

(deg) 

(deg) 

Lat. 

(deg)  Long. 

Computed  Dip 

West 

East 

East 

(deg) 

Adak            AD 

51.90  N 

176.60 

183.40 

47.22  N 

239.99 

63.1  N 

Bogota          BG 

04.50  N 

74.20 

285.80 

15.95  N 

354.61 

23.8  N 

Boulder          BC 

40.00  N 

105.30 

254.70 

48.85  N 

316.44 

68.3  N 

Conception       CP 

36.60  S 

73.00 

287.00 

25.12  S 

356.15 

35.5  S 

Grand  Bahama      GB 

26.60  N 

78.20 

281.80 

37.93  N 

349.56 

59.8  N 

Huancayo         HU 

12.05  S 

75.30 

284.70 

00.62  S 

353.81 

1.3  N 

Jamaica         JA 

18.00  N 

76.80 

283.20 

29.38  N 

351.48 

50.4  N 

Manila          MiN 

14.70  N 

121.10 

03.37  N 

189.78 

14.5  N 

Maui            MA 

20.80  N 

156.50 

203.50 

20.86  N 

268.11 

38.8  N 

Mexico  City       MX 

19.40  N 

99.70 

260.30 

29.14  N 

327.08 

46.4  N 

Okinawa          OK 

26.30  n 

127.80 

15.25  N 

195.58 

36.6  N 

Point  Arguello    PA 

35.60  N 

120.60 

239.40 

42.15  N 

300.78 

61.0  N 

Talara          TA 

04.60  S 

81.30 

278.70 

06.64  N 

347.66 

13.1  N 

Wallops  Island    WP 

37.90  N 

75.50 

284.50 

49.31  N 

352.12 

69.9  N 

White  Sands      WS 

32.30  N 

106.50 

253.50 

41.10  N 

316.94 

60.5  N 

-These  va 1 ues  will 

not  agree  wi 

th  those 

of  Report 

UAG-38, 

Master  Statior 

List  for 
December 

Solar-Terrestrial 

Physios  Data 

at  WDC-A  for  Solar-Terrestrial  Physios. 

\31h.      Different 

field  models 

were  uset 

i. 

Height  Profile   of 

Electron  Dens' 

Table   3.      Sample 

printout;   Single   True 

'.ties 

BOULDER,    COLO.  73    06    15    0700 

FOR    SAMPLE         F=0 . 560000  FH=1. 429120 

LOGFN  H*  H  Z* 


SIN    1=0.924546 

Z**  INTEGRAL 


9.9  0  50  EXTRAP 
0.005    90  O 


0.015 
0.250 
0.350 
0.370 
0.400 
0.420 
0.430 
0.440 
0.450 
0.470 
0.500 
0.530 
0.560 
0.580 
0.593 
0.603 
LOG  FM=0.6067 


90  O 

107  O 

110  O 

119  O 

120  O 
131  O 
150  O 


90.00 

90.00 

90.00 

94.51 

98.78 

99.99 

102.04 

132.70 

105.37 


218  OV  109.51 
190  O   112.52 


181  O 

190  O 

210  O 

251  O 

307  O 

395  O 

586  O 


117.49 
124.27 
132.26 
14  3.16 
154.09 
165.46 
181.75 


0.000  0.00 

0.000  0.00 

0.000  0.00 

20.915  44.50 

17.801  -000015.57 
42.475  616.85 
25.630  -000280.75 

57.802  804.29 
108.706  2545.21 
305.188  9824.13 
141.785  -000882.00 
106.505  -000882.00 


HMAX=199. 86 


119.715 
146.497 
217.021 
329.392 
545.043 
1084.418 
S=    69.527 


220.17 

446.36 

1175.41 

2809.28 

8294.26 

26968.74 

SHMAX= 


0.000 

0.000 

0.000 

1.386 

3.302 

4.092 

5.576 

6.952 

8.417 

12.236 

15.134 

20.247 

28.117 

38.761 

55.486 

74.257 

95.306 

127.119 

163.608 


EL/Cu 

00805+003 
26919+004 
32901+004 
92217+004 
21623+004 
81596+004 
82576+004 
8.58078+004 
8.98518+004 
9.40864+004 
,85205+004 
08026+005 
24030+005 
42405+005 
63503+005 
79278+005 
90339+005 
99309+005 
.02697+005 


Line    1    —    Station,    Year,    Month,    Day,    Hour,      Minute    (local    standard   meridian    time    - 

Line   2    --    Project    designator    (for   sample),    Geomagnetic    field   strengths    at   ground 
level     in    gauss    (F) ,    gy rof requency    at   200   km    (2.55xF)     in   Megahertz    (FH) , 
the   sign   of    the    local    geomagnetic   dip   angle    (SIN    I). 

Line   3   --    (gives    the    column    head i ngs) --1 ogar i thm  of    the    plasma-frequency    (LOG    FN), 
virtual    height    (H-)  ,    true   height    (H) ,    slope   of    the   profile    (Z-) ,    change 
in    the   slope    (Z**) ,    integrated   electron    density    (N10    electrons/cm2) 
(INTEGRAL),    and    the   electron    density  with    the  exponents    as    shown    (EL/CC) 

Last   line   --    Logarithm   of    the   maximum   plasma    frequency    (LOG    FM) ,    height   of    the 

maximum   plasma    frequency    (HMAX) ,    the   quarter   thickness    (S   or   sometimes 
SCAT),    the    integrated   electron   density    up    to   HMAX    (SHMAX) ,    the   electron 

density   at   HMAX. 
"V"    in   H-    column    represents    a    critical    frequency   and   "0"    indicates    the   ordinary 
trace  was    used    to  obtain    the    virtual    heights. 
All    heights    in    kilometers. 


LMT) 


Table  4.      Sample  Printout;   Card  Images  of  Five  Interpolated  Electron  Density  Profiles 


540000 
540000 
540000 
500000 
500000 
500000 
500000 

5 JOQ64 

5J0064 
5J0064 
520132 
523132 
^20132 


BC640316 
BC640316 
BC640316 
BC640316 
BC640316 
8C64Q316 
BC640316 
BC640316 
BC640316 
BC640316 
BC640316 
8C64Q316 
BC640316 


1700100 
1700110 
1700120 
1715110 
1715120 
1730110 
1730120 
1745100 
1745110 
1745120 
1800100 
1800110 
1800120 


0000 
0  362 
2790 
0312 
2830 
0270 
2490 
0000 
0236 
3370 
0000 
0197 
2650 


00000 
00511 
03690 
00459 
03750 
00389 
03510 
00000 
00371 
03960 
00000 
00246 
03370 


000000000 
005500059 
043800475 
005070054 
044200476 
004230045 
043300469 
0000  0  0000 
004C20046 
042600434 
000030000 
002770031 
037830388 


00000 
10066 
00484 
40057 
00484 
00048 
00473 
00000 
90055 
00434 
00000 
60037 
00388 


00000000 
60078000 
00485000 
90063200 
00485000 
50053200 
00474000 
00000000 
6C068000 
00000000 
00000000 
00045400 
00000000 


0000000 
9230112 
0000000 
7950103 
0300000 
6010071 
0000000 
0000000 
8760121 
0000000 
0000000 
6010085 
0000000 


000000 
001400 
000000 
001340 
000000 
80C971 
000000 
000000 
001760 
0C0000 
000000 
601270 
000000 


00079 
01870 
00000 
019  30 
00000 
01530 
00000 
03055 
02570 
00000 
00049 
01870 
00300 


625  0325 
00OOOO0OGO 
090239C26M 

625  0329 
090239025Q 

618  0291 
090236022J 

591  0313 
0000000000 
090228022L 

559  0274 
3000000300 
090229017Q 


Columns    1    througn    6    --    give    the    internal    job    identification. 

Columns    7    and    8    --    give    the    two-letter   station    code    (see    Table    2) 

Columns    9    through    18    —    give    the   year,    month,    day,    hour,    minute     (LMT). 

Column    19    --   gives    the   profile   quality    grade,    see   section    5.1. 

Column    20    --    is    a    key    digit   which    identifies    the    height    range    for    the    data    in    columns    21 

through    70    (0=0    to    90    km,     1    =    100    to    190    km,    2    =    200    to    290    km,    etc.). 
Columns    21    through    70    are    divided    into    ten    five-column    sections.       Each    gives    the    electron 

density    in    10      electrons/cc    in    increments   of    ten   km. 
Columns    71    through    80    --    of    the    first    card    of   each    profile    give    the    foF2    of    ionogram    in 

columns    71"73    and    the    quarter    thickness    (SCAT)    of    the    F2    layer    in    columns 

7^-80. 
Columns    71     through   80    --    on    the    last    card    of    each    profile    q i ve    the    minimum    height    of    the 

profile    (HMIN,     in    columns    71     through    73),    the    maximum   height    of    the    profile 

(HMAX,     in    columns    Jh    through    76)    and    the    integrated   electron    density    up    to 

the    height    of   maximum   plasma    frequency    (SHMAX,     in    columns    77    through    80) . 

The    number    in    column    80   of    the    last    card    receives    an    overpunch    to    flag    it 

as    being    the    final    card   of    the   profile. 
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Figure  1.      Electron  Density  Profile  -from  Ionogram  recorded  at  Boulder  on  December  10,    1964, 
at  midnight    (BC  641210   0). 


Table  5.      Interpolated  Electron  Density  Profile  Output 
(Electron  Densities  in  10z  EL/CC) 
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3.        THE  LERFALD  ANALYSIS   SYSTEM 
3.1     Film  Review 

Ionograms  on  35mm  film  for  the  15  stations  of  Table  2  for  each  of  the  4  seasons  in  1964  (Table  1) 
were  retrieved  from  the  archives  of  World  Data  Center  A  for  Solar-Terrestrial  Physics  and  sent  to  the 
EDS  Ionogram  Quality  Review  Group  for  comment.  Each  reel  was  scanned  to  note  information  to  aid  inter- 
pretation and  digitization.  At  least  one  sample  ionogram  print  was  made  for  each  station.  This  print 
was  annotated  to  reveal  fiducial  points,  distortions  due  to  equipment  misalignment  and  possible  diffi- 
<lul^r^sjrin  interpretation.  All  of  these  specific  comments  and  guidelines  by  EDS  specialists  were  sent 
to  SEL  for  use  by  its   contractor. 

Generally,  only  stations  with  reputations  for  having  well-maintained  equipment,  easy-to-scale  35mm 
film  ionograms  and  midlatitude  geographic  locations  were  selected  for  this  project.  The  stations  known 
to  be  easier  to  handle  were  sequenced  to  the  digitizer  ahead  of  the  more  troublesome  ones 


3.2  Measuring  the  Data 

In  order  to  automatically  process  the  ionogram  data,  sufficient  information  must  be  extracted  from 
each  ionogram  to  define  the  h'(f)  traces  and  to  transform  them  into  digital  format.  An  automatic  trace- 
following  laser  beam  measuring  machine  (called  the  SWEEPNIK)  [Lerfald  et  al .  ,  1974]  was  interfaced  with 
an  online  computer.  An  operator  with  manual  override  capabilities  supervised  the  digitizing  to  insure 
the  correct  trace  was  beinq  followed.  This  system  utilized  the  advantages  of  the  inherent  accuracy  of 
the  laser,  the  speed  of  the  computer,  and  the  pattern  recognition  capabilities  of  the  operator. 

The  resulting  digital  magnetic  tape  contained  station  and  frame  identification,  fiducial  points,  x-y 
coordinate  values  of  the  ionogram  trace  (in  cm)  and  coded  comments.  A  small  degree  of  smoothing  was  pro- 
gramed into  the  laser  optics  to  prevent  undesirable  jitter  in  the  digital  data. 

3.3  Computer  Reconstructing  and  Formating  of  the  Data 

The  measured  data  were  then  ready  to  proceed  throuqh  a  series  of  subroutines  to  prepare  them  for 
processing  by  the  HOWE  inversion  program  [Howe  and  McKinnis,  1967].  The  subroutines  are  briefly 
described  as  fol lows : 

1)  SURVEY  -  lists  station  code,  date,  time,  number  of  data  points,  number  of  fiducial  points, 

the  first  and  last  frequency,  tne  difference  in  the  first  and  last  frequency,  the  difference  between  the 
frame  fiducial  and  the  master,  and  operator's  comments  for  each  frame. 

2)  START  -  assigns  the  D-region  model  to  be  used  in  the  profile  calculation. 

3)  ELAYER  -  assigns  E-layer  data  points  from  a  model  when  none  is  scaled. 

4)  X-YPICK  -  selects  the  digital  data  points  to  be  used  in  the  final  profile  calculation. 

5)  VPOIHT  -  assigns  "V"  designations  to  discontinuities  at  layer  peaks. 

6)  TRANSFRM  -  converts  the  x-y  coordinate  values  to  h'(f)  values. 

7)  FIVECHAR  -  checks  tne  one-hop  to  two-hop  ratio. 

8)  TYPE2  -  formats  data  for  HOWE  input. 

3.4  Constructing  the  "Unseen"  Ionization  Model 

The  "unseen"  ionization  must  be  taken  into  account  to  correct  the  virtual  height  of  the  first 
observed  point  on  the  ionogram  to  a  "true"  height.  It  is  represented  by  a  model  which  is  applied  to 
the  observed  ionogram  data  in  two  steps.  The  first  step  is  to  establish  the  electron  density  at  90  km 
(below  which  echoes  are  seldom  observed).  Next,  virtual  heights  points  must  be  selected  to  connect 
the  first  observed  data  with  the  90  km  points. 

3.4.1.  F(90)  Model 

The  F(90)  Model  of  the  electron  density  at  90  km  is  due  to  Smith  [1974]  and  in  effect  takes 
account  of  the  diurnal  variation  of  the  D-region  ionization  by  the  simple  relationship  involving  the 
cosine  of  the  solar  zenith  angle  by  day  and  linear  relaxation  in  the  late  afternoon  and  evening  (see 
Figure  2).  This  model  was  designed  to  match  radio  propagation  experience. 


hours  before  sunrise  or 
after  sunset 


cosine  of  the 
solar  zenith  angle 


Figure   2.      Smith  F(90)   Model  of  Diurnal   Variation  of  D  Region. 


Lerfald's  attempt  at  expressing  this  curve  in  computer  code  is  explained  in  SEL-34  [Jurgens  et  al. 
1974].  His  model  was  limited  to  latitudes  lower  than  approximately  55°,  thus  to  those  stations  having' 
at  least  5  hours  of  darkness.   However  the  model    has   been  extended  recently  for  global   coverage   [Smith, 

3.4.2.       fie  E-  and   F-region  Model 

he  E-  and   F-region  Model    begins  with  a   predicted  value  for  the  E-layer  critical    frequency   (foE) 
as       function  of  time  of  day,   season,   sunspot  cycle  and  geographic  location.     This   value  is   taken  from 
tr      empirical    E-layer  climatology  used  for  ionospheric   radio  propagation  predictions   [Roberts  and 
R    ,ich,   1971].     The  foE  is  zenith  angle  dependent  by  day.     Nighttime  values  are  from  limited  observa- 
t    ins  and  radio  propagation  experience.     The  true  height  profile  of  the  E  layer  is  assumed  to  be 
i    jghly  parabolic,   and   virtual    height  points   for  the  model   were  selected  to   reflect  this. 

The  recording  on  the  ionogram  of  the  virtual  height  of  critical  frequency  is  dependent  upon  trans- 
mitter power  and  receiver  sensitivity.  In  the  vicinity  of  foE,  for  example,  it  is  common  for  the  trace 
to  be  so  faint  that  neither  foE  nor  its  virtual  height  can  be  actually  measured,  yet  a  numerical  height 
value  must  be  provided  to  the   "HOWE"   program. 

For  modeling  an  unseen  E  layer,   Jurgens  et  al.    [1974]    usually  set  the  virtual    height  of  the  E-layer 
cusp  (foE)  at  50  km  above  the  first  scaled  height  in  the  F  layer.     An  exception  to  this  practice  was 
when  the  lower  end  of  the  F  layer  trace  had  a  zero  or  positive  slope.     In  this  case  the  virtual    height 
of  foE  was  set  equal    to  the  height  of  the  first  scaled  point  on  the  F  layer.     Serious   difficulties  occur 
in  the  calculation  of  the  electron  density  distribution  when  the  foE  virtual    height  is  close  to  the 
virtual    height  of  the  first  F-layer  point.     Here  is   the  only  major  difference  between  Jurgen's  work  and 
the  model   adopted  by  EDS   (Figure  3).     Because  of  the  built-in  continuity  equations   of  the  "HOWE"   program, 
it  makes   little  or  no  difference  whether  the  virtual    height  is  excessively  high  or  not.     EDS  arbitrarily 
assigns  400  km  as  the  virtual   height  of  the  E-layer  cusp  in  an  effort  to  make  the  profile  computation 
more  consistent. 
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Figure  S.      E-  and  F-Region  Model  used  by  EDS. 

Figure  3  shows   the  E-   and  lower  F-region  EDS  model   with  a  400  km  f oE .      It  should  be  noted  that  all 
E-layer  heights   shown  are  adopted  representative  virtual    heights   and  distributed   in   frequency  relative 
to  the  foE.     The  shape  of  the  E  layer  is  specified  by  three  height-frequency  points,  as  follows: 


(118  km  at  0.86  MHz  x  foE),    (124  km  at  0.93 
(150   km  at  0.98  MHz  x  foE)  . 


\z   x  foE) ,  and 


The  shape  of  the  model  E  layer  is  best  suited  to  low  sunspot  years,  midlatitude  stations  and  magnetically 
quiet  summer  midday  ionograms.  No  allowance  has  been  made  for  the  height  variation  caused  by  temporal 
or  magnetic  changes.  The  heights  are  not  adjusted  other  than  placing  them  at  86%,  93%,  and  98%  of  the 
scaled  or  assumed  foE.  The  foE  frequency  was  varied  in  relation  to  sunspot  number,  season,  latitude, 
and  time  of  day. 

The  E  region  is  joined  to  the  measured  F-region  data  by  a  fi rst-order  approximation  to  an  exponen- 
tial curve. 


3.5  The  Inversion  Method 

The  problem  faced  here  is,  given  the  virtual  height  as  a  function  of  frequency  obtained  by  the 
ionosonde  technique,  to  calculate  the  electron  density  distribution  as  a  function  of  real  height.  This 
is  a  troublesome  problem  which  can  never  have  an  unambiguous  answer,  because  the  inversion  integral  to 
convert  h'(f)  into  h(f)  requires  complete  information  about  h'(f)  over  all  frequencies,  and  this  cannot 
be  obtained  by  the  ionosonde  technique.  An  infinite  number  of  electron  density  profiles  can  produce 
the  same  incomplete  ionogram.  Thus,  certain  assumptions  must  be  made  about  the  nature  of  the  electron 
density  distribution  before  a  calculation  can  be  made  with  the  available  h'(f)  information.  To  begin 
with,  some  assumption  must  be  made  about  the  effect  of  underlying  ionization  with  densities  less  than  that 
which  can  reflect  the  lowest  frequency  seen  in  the  echoes  on  the  ionogram.  Another  common  assumption 
is  that  the  electron  density  increases  monotonically  with  height.  A  detailed  description  of  the  step- 
by-step  procedure  of  the  calculation  is  given  by  Wright  [1967]. 

In  the  method  used  to  calculate  the  data  base  presently  under  discussion  it  was  assumed  that  dN/dh 
is  greater  than  0  until  foF2  is  reached,  at  which  time  dN/dh  becomes  0.  This  monotonic  condition  can 
sometimes  be  proved  by  inspection  of  the  ionogram.  However,  more  often  it  is  not  confirmed  and  the 
existence  of  a  "valley"  between  the  E  and  F  region  is  a  very  real  possibility  even  though  conventional 
ionograms  give  no  direct  measurement  of  the  morphology  of  that  region.  The  HOWE  program  does  allow  for 
at  least  partial  correction  of  this  valley  situation.  Such  a  correction,  however,  was  not  included  in 
the  Lerfald  automatic  processing  adaptation. 

It  is  also  impossible  to  measure  from  an  ionogram  the  exact  value  of  the  foF2  because  the  virtual 
height  is  infinite  at  that  frequency.  Customarily  foF2  is  scaled  by  estimating  the  asymptotic  position. 
The  HOWE  program  does  not  use  the  scaled  value  of  foF2  but  instead  uses  a  systematic  method  of  extrapo- 
lating the  observed  trace  to  foF2  and  of  extrapolating  to  the  height  at  which  the  maximum  electron 
density  occurs.  This  is  also  described  by  Howe  and  McKinnis  [1967]. 

3.6  Internal  Program  Rejection  Parameters 

An  important  feature  incorporated  into  the  HOWE  program  is  the  ability  to  flag  bad  profiles.   In 
addition  to  normal  computer  input  data  checks,  this  program  also  checks  the  output  data  for  profile 
parameter  variations  beyond  a  subjectively  established  reasonable  limit.  The  parameters  checked  are 
the  profile  slope,  integrated  electron  density,  height  of  maximum  electron  density,  height  of  minimum 
electron  density,  the  quarter  thickness  of  the  parabolic  nose,  the  parabolic  fit  of  the  peak,  the 
magnetic  field,  and  the  standard  deviation  of  the  "start"  extrapolation.  These  parameters  are  also 
listed  on  the  profile  output  so  a  more  detailed  check  can  be  made.  In  the  hand-processing  of  profile 
data  this  detailed  checking  is  an  important  step  in  the  quality  control.  It  is  also  one  of  the  most 
difficult  things  for  the  inexperienced  person  or  for  a  computer  to  do.   In  the  Lerfald  SEL  automated 
system  this  step  was  neglected  altogether. 

4.   APPLICATION  OF  THE  LERFALD  UNDERLYING  IONIZATION  MODEL 

To  demonstrate  the  model  application  techniques,  a  typical  afternoon  ionogram  is  shown  in  Figure  4 
in  several  degrees  of  observed  detail,  and  the  model  is  applied  appropriately  to  each  case. 

Figure  4a  is  the  easiest  ionogram  to  process.  Heights  are  observed  and  can  be  measured  down  to 
about  100  km  at  fmin,  so  this  ionogram  is  nearly  complete,  since  there  is  little  ionization  below 
100  km.  The  Lerfald  model  needs  only  to  supply  the  90  km  point  (as  determined  by  Smith  [1974]}.  For 
monotonic  cases  the  profile  resulting  from  the  calculation  of  these  virtual  height  data  is  as  accurate 
as  is  possible  from  the  method  used. 

All  of  the  F  region  and  a  definite  indication  of  foE  can  be  seen  in  Figure  4b.  In  this  case  three 
height-frequency  values  for  the  E  layer  (their  frequency  depending  on  the  value  of  foE  adopted  from  the 
indication  on  the  F-region  trace)  and  one  for  the  D  region  are  supplied  by  the  model.  The  accuracy  of 
the  resulting  profile  depends  on  how  well  these  three  points  describe  the  real  shape  of  the  "unseen" 
E  region  and  the  D-region  electron  distribution  (as  described  above). 

Figure  4c  shows  an  ionogram  which  gives  an  indication  of  underlying  ionization  by  the  minor  cusp  at 
the  lower  end  of  the  F  region,  but  little  else.  In  this  case  an  foE  is  calculated  for  the  season,  hour 
and  sunspot  number  [Roberts  and  Rosich,  1971]  and  set  at  a  virtual  height  of  50  km  above  the  first  scaled 
point  in  the  F  region.  The  resulting  profile  has  additional  inaccuracies  caused  by  prediction  and  model 
connection  errors. 

Progressing  a  bit  further  with  the  absence  of  necessary  information,  as  shown  in  Figure  4d,  we  note 
that  none  of  the  retardation  caused  by  underlying  ionization  is  observed  on  the  ionogram.  Lerfald  de- 
cided to  handle  this  situation  by  assigning  the  virtual  height  of  foE  to  be  the  same  as  the  height  of 
the  first  scaled  point  of  the  F  layer.  As  described  earlier  (see  Sec.  3.4.2),  this  can  lead  to  trouble 
in  the  profile  calculation  by  greatly  lowering  the  computed  heights.  A  better  approach  would  be  to 
assign  an  arbitrary  height  of  about  400  km,  as  in  the  EDS  modeling  technique. 
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Figure  4.      Sample  Afternoon  Ionogram  with  Successive  Deletion  of 
Lower  Frequency  Echoes  to  Illustrate  Model  Application 
as  Observed  Detail  Decreases .      0-trace  Only  is  Shown. 


Finally,  the  absence  of  lower  level   data  (Figure  4e)   is  such  that   it  is  doubtful   that  the 
information  on  the  ionogram  warrants  an  attempt  to  compute  a  profile.     Generally  if  the  slope 
(dh'/dlog  (f^))  is  much  greater  than  50,  the  results  of  an  attempt  to  calculate  a  profile  by  this 
method  are  doubtful . 

These  examples  illustrate  a  range  of  conditions  experienced  in  practice  in  determining  N(h)  pro- 
files from  actual   ionograms,   in  this  case  daytime  ionograms.     The  profiles  calculated  in  these  five 
examples  are  tabulated  in  Table  6.     There  is  fair  consistency  among  the  profiles  calculated  from  (a): 
(b)  and  (c)  of  Figure  4.     The  apparent  degradation  when  so  much  information  is  missing  can  be  seen 
in  the  results  from  (d)  and  especially  (e). 

Table  6.      Profile  Calculations  of  Ionograms  in  Figure  4 
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5.      PROFILE  QUALITY  ASSESSMENT 

Assessing  the  quality  of  an  electron  density  profile  is  not  a  straightforward  task  because  of  the 
da  oinus°n  he  1p  hoHriab  "^  IS6  ^aliVePends  upon  the  accuracy  and  completeness  of  tie     onogram 
vals  f      m?ft  C2S  US       a       the  C3re  ta,ken  in  the  comPutatTon.     Since  there  were  short  time  inter- 
vals  (15  mm)  between  Tonograms  processed  for  this  data  base  and  since  changes  in  the  ionosphere  were 
generally  smooth  over  such  time  scales,  we  decided  the  most  meaningful   test  wou  d  be  to  check  the 
n,nf  Cy  °f  the  Proflles  from  one  P^file  to  the  next.     This  will   not  necessarily  be  a  test  of 

Ki'S'SSiitaj  ?hetph?of?li:'StenCy  "'^  Wh1Ch  the  Pr°JeCt  eX6CUted  the  -thodVanalvl^which 

bv  rhIrkinnSiLWHle-^ed  t0/heck  t°r  consistency.     Each  profile  in  the  entire  data  base  was  graded 
calculatld9f^  Si         ?n!°f  a  Pr0fl1!  Parameter  from  the  average  value  of  the  profile  parameter 
the  t  » of  Z  Si  tEI°fllf-iPTeKSed  inM  ,that  m°nth  f0r  that  Station  and  month  within  one  hour  of 
HnL     w     I    fy  °f  the>  Pr0flle  t0  be  graded-     In  addition  to  grading  the  entire  data  in  this  way,  a 
nrlrln       '^ ,T  ""^  °f  pr0files  fbr  selected  times"     The  latter  were  hand-processed  by  the  same 
procedures  used  to  compute  the  profiles  in  the  automated  mode.     These  are  discussed  in  the  following 
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5.1  Grading  the  Profiles 

Each  profile  was  assigned  a  grade  of  acceptable  (grade  1),  marginal  (grade  2),  unacceptable  (grade 
3)  and  rejected  (grade  9).  The  scheme  involves  combining  individual  grades  based  on  deviation  of  each 
of  several  profile  parameters  from  their  average  in  stated  ranges.  Table  7  gives  the  details  and  the 
relative  weights  given  to  the  grading  by  the  various  parameters. 

Table   7.      Variation  and  Weights  of  Profile  Parameters 


Profile      Grade  1     Grade  2      Grade  3 

Parameter    acceptable    marginal     unacceptable    weight 


foF2 

0  to  14% 

15  to  29% 

30  &  above 

2 

Hmax 

0  to  9% 

10  to  19% 

20  &  above 

3 

Scat 

0  to  49% 

50  to  99% 

100  &  above 

1 

Shmax 

0  to  24% 

25  to  49% 

50  &  above 

2 

N(200) 

0  to  19% 

20  to  39% 

40  &  above 

2 

The  variations  are  deviations  with  respect  to  the  corresponding 
two-hour  average  for  the  parameter  concerned. 

Profiles  rejected  altogether  and  released  from  the  data  base  were 
assigned  a  numerical  grade  "9". 

Spot  checks  were  made  of  profiles  which  were  graded  two  or  three  in  order  to  find  the  reason  for 
the  inconsistencies.  In  most  cases  this  search  revealed  that  either  there  had  been  an  incorrect  inter- 
pretation in  the  ionogram  or  the  incorrect  edge  (following  rather  than  leading)  of  the  echoes  had  been 
digitized.  Interpretation  of  an  ionogram  is  relatively  easy  during  magnetically  quiet  periods  when  the 
echoes  are  being  reflected  from  a  nearly  overhead  horizontally  stratified  ionosphere.  However,  this  is 
not  the  case  most  of  the  time,  and  there  are  times  when  even  experts  argue  over  the  proper  interpre- 
tation. 

Figure  5  shows  an  example  of  an  ionogram  from  which  a  poor  profile  was  produced,  A  good  profile 
would  have  been  possible  had  the  interpreter  been  somewhat  more  experienced.  The  first  reason  the 
ionogram  resulted  in  an  unacceptable  profile  was  that  the  second  reflection  of  the  c-type  sporadic  E 
layer  was  digitized  as  echoes  from  the  F  region.  The  second  reason  the  ionogram  resulted  in  a  bad 
calculation  was  incomplete  digitizing  of  the  data.  Only  the  E  and  Fl  layers  were  digitized;  the 
distinct  F-region  trace  at  a  great  height  was  overlooked.  Seldom  would  an  experienced  scaler  have  made 
these  errors.  In  this  example  one  can  also  see  the  leading  edge  problem  described  earlier. 

Each  profile  in  the  data  base  has  been  assigned  a  grade,  and  this  appears  along  with  the  profile 
stored  on  magnetic  tape.  It  is  also  indicated  in  the  profile  index  of  which  Table  8  is  a  sample.  The 
index  for  all  profiles  in  the  data  base  appears  in  Conkright  [1977]. 

5.2  Comparison  of  Hand-Analysed  Versus  Mass-Produced  Profiles 

A  detailed  analysis  of  profiles  for  selected  times  has  been  conducted  to  compare  the  results  of  the 
automated  system  with  those  obtained  by  hand  processing  the  ionograms.  Both  procedures  followed  essen- 
tially the  same  rules  regarding  model  application  and  the  inversion  calculations. 

The  exception  was  in  the  way  the  model  for  underlying  ionization  was  applied  to  the  (h'f)  point  at 
foE  as  discussed  above.  In  the  hand  processing  the  EDS  scheme  of  putting  h ' =400  at  foE  was  used   Also 
there  proved  to  be  some  small  differences  in  the  climatol ogical  prediction  used  for  the  value  of  foE. 
In  most  cases  these  exceptions  made  no  significant  difference  for  values  obove  the  minimum  observed 
frequency  of  the  normal  layer  (Fl  or  F2).  Data  below  the  lowest  frequency  observed  from  the  normal 
layer  on  an  ionogram  should  be  disregarded  in  any  case,  since  they  are  only  model  data. 

The  ranges  for  the  different  grades  were  rather  arbitrarily  determined  by  study  of  the 
ionograms  discussed  in  this  section.  The  profiles  from  the  hand  analysis  were  considered  to  be  the 
"truth".  Parameters  from  these  were  compared  to  the  corresponding  parameter  from  the  mass-produced 
profiles  in  each  of  several  ranges  of  the  deviation  from  average.  The  agreement  between  hand  and  mass 
results  was  reasonably  good  in  the  ranges  like  adopted  for  grade  1,  and  deteriorated  at  larger  ranges. 

Similarly  the  weights  assigned  to  each  parameter  for  forming  the  grade  level  for  the  profile  cal- 
culation as  a  whole  are  somewhat  arbitrary.  Discrepancies  in  "Hmax"  between  hand-and  mass-produced 
profiles  seemed  to  be  a  rather  better  or  more  sensitive  discriminator  of  consistent  profiles  while 
"Scat"  was  less  sensitive  than  the  others. 
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A  summary  of  the  grades   assigned  to  the  18,226  profiles   in  the  data  base  is  shown   in  Table  9. 

The  station  and  time  selections  for  these  analyses  were  based  on  four  easy-to-handle  stations 
with  a  history  of  good  performance   (Point  Arguello,  White  Sands,  Manila,   and  Boulder)   and  on  periods 
of  appreciable  contrasts,   that  is,   noon,  midnight,   sunrise   (0400,0500  and  0600  local   time)   and 
sunset   (1600,   1700,    1800  local    time). 
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Figure   5.      Boulder  0920 (LT J    March   17,    1964.      Ionogram. 
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Table  9. 

Summary 

of  Profile 

Grade 

s 

Station 

1964 
Month 

if  o 

1 

f  Pro 

2 

files  graded: 
3      9 

Station 

1964 
Month 

if   o 
l 

f  Prof 
2 

iles  gi 
3 

aded : 
9 

Adak 

May 
Jun 
Sep 
Dec 

247 
201 
502 
578 

11 

9 

21 

51 

1 
0 

0 
1 

1 
9 
1 
0 

Maui 

Mar 
Jun 
Sep 
Dec 

210 
177 
344 
389 

30 
28 
70 
95 

0 
0 
0 
0 

7 

6 

15 

16 

Bogota 

May 
Jun 
Sep  ■ 
Dec 

229 
216 
366 
419 

31 
27 

60 
37 

0 
0 

0 
0 

0 
0 
0 
0 

Mexico  City 

May 
Jun 
Sep 
Dec 

59 
160 
316 
598 

0 
21 
38 
37 

0 

1 
0 

1 

0 

4 
1 
2 

Boulder 

May 
Jun 
Sep 
Dec 

275 
168 
406 
641 

1 
9 

46 
39 

0 

0 
0 
0 

1 
5 
7 
0 

Okinawa 

May 
Jun 
Sep 
Dec 

201 

37 

396 

504 

20 

4 

56 

92 

0 
0 
2 
1 

0 
0 
0 
2 

Concepcion 

May 
Jun 
Sep 
Dec 

195 

215 
371 
439 

6 

10 
35 
31 

0 
0 
0 
2 

0 
0 
0 
3 

Pt.  Arguello 

May 
Jun 
Sep 
Dec 

117 
147 
208 
117 

5 

18 

8 

198 

0 
2 
0 

1 

3 
7 
6 
8 

Grand  Bahama 

May 
Jun 
Sep 
Dec 

227 
145 
343 
621 

14 
12 
16 
52 

0 
0 
0 
0 

0 
1 
1 
2 

Talara 

May 
Jun 
Sep 
Dec 

182 
186 
304 

476 

29 

29 
60 
51 

0 
0 

1 

0 

0 
0 
0 
0 

Huancayo 

May 
Jun 
Sep 
Dec 

0 
163 
196 
469 

0 
41 
26 

40 

0 
0 
0 

1 

0 
2 
5 
7 

Wallops 

Is. 

May 
Jun 
Sep 
Dec 

0 

0 
0 
0 

0 
0 
0 
0 

0 
0 
0 
0 

0 
0 
0 
0 

Jamaica 

May 
Jun 
Sep 
Dec 

53 
138 
358 
324 

0 
16 
29 
19 

0 
0 
2 
0 

0 
6 

1  1 
6 

White  Sands 

May 
Jun 
Sep 
Dec 

157 
163 

416 
600 

7 

5 

61 

42 

0 
0 
0 
0 

0 
2 
2 
6 

Manila 

May 
Jun 
Sep 
Dec 

166 
155 
358 

477 

34 

30 

66 

103 

0 
0 

1 
1 

0 
2 

1 
0 
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In  Figure  6a,  b,  c,  and  d  can  be  seen  the  graphic  results  of  the  Boulder  portion  of  the  data 
comparisons.  Studying  the  comparison  of  profiles  for  the  selected  stations  and  times  reveals  several 
generalities:  1.  When  the  scaled  ionogram  values  are  the  same  by  both  the  hand  and  automated  methods, 
and  the  foE  can  be  seen,  the  resulting  profiles  are  essentially  identical;  2.  When  data  cannot  be 
scaled  to  frequencies  at  or  below  foE,  errors  attributed  to  the  model  predictions  and  application  can 
result.  Differences  in  interpretation  of  the  ionograms  varied  from  small  defensible  differences  of 
opinion  to  gross  mistakes  by  the  interpreter. 

Table  10  presents  a  general  idea  of  the  statistical  differences  between  the  hand  and  automated 
methods  of  processing  the  same  ionograms.  The  statistical  comparison  is  divided  into  noon,  midnight, 
sunrise,  and  sunset  periods.  The  rationale  behind  this  type  of  subdivision  was  that  the  ionogram  data, 
regardless  of  season  or  station,  would  be  more  or  less  homogeneous  within  each  of  these  divisions.   In 
other  words,  the  noon  ionograms  are   the  most  likely  to  show  the  largest  amount  of  information  useful  in 
the  N(h)  inversion,  (E,  Fl  and  F2  regions);  and  the  midnight  ionograms  will  most  likely  show  the  least 
(F  region  only)  information.  Also,  the  monctonic  assumption  will  be  more  comparable  within  these 
divisions. 

Table   10.      Comparison  of  Profile  Parameters 


AHmax 


ANmax 


AScat 


AShmax 


AN 


200 


BOULDER 
SUNRISE 

NOON 
SUNSET 
MIDNIGHT 


3.0 
2.4 
0.0 
2.1 


10.7 
0.8 

11.1 
6.3 


4.1 

11.8 

2.1 

1.0 


12.9 

12.7 

7.3 

2.1 


6.0 

7.9 

10.0 

17.2 


PT  ARGUELLO 
SUNRISE 

NOON 
SUNSET 
MIDNIGHT 


2.6 

NO  DATA 

12.0 

6.9 


1.5 

NO  DATA 

22.1 

3.2 


4.5 

NO  DATA 

9.5 

15.9 


9.8 

NO  DATA 

10.7 

15.7 


23.3 

NO  DATA 
22.0 
38.5 


WHITE  SANDS 
SUNRISE 

NOON 
SUNSET 
MIDNIGHT 


7.6 
1.1 
1.2 
1.7 


2.1 

1.1 
0.2 
2.0 


12.9 
25.6 
24.5 


1.9 
7.2 

5.3 
10.0 


3.1 

0.4 

6.9 

15.4 


MANILA 
SUNRISE 
NOON 
SUNSET 
MIDNIGHT 


13.9 
9.7 
2.3 
1.6 


2.3 
7.1 
0.4 
0.1 


9.7 
25.6 
11.5 
11.9 


1.1 

20.2 

4.5 

0.0 


1.6 

2.8 
10.8 
16.8 


Average  differences  between  profile  produced  by  hand  (EDS)  and 
semiautomatic  (SEL)  methods.   All  values  in  %. 


Time  and  funding  permitting,  large  portions  of  trie  Lerfald  Data  Base  should  be  redone  to  correct 
for  the  gross  deficiencies  that  undoubtedly  biased  an  analysis  of  this  kind  (see  Figure  6c,  Boulder 
September  data  as  an  example).  A  consistency  comparison  such  as  we  used  cannot  flag  a  constant  error  no 
matter  how  large.  From  a  cursory  survey  of  the  hand  vs.  automated  analysis,  it  would  seem  tiiat  all  or 
the  September  data  from  Boulder  had  an  error  of  this  type  and  therefore  should  either  be  rerun  or 
disregarded. 
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5.3  Conclusions  of  Profile  Comparisons 

The  conclusion  one  can  draw  from  this  quality  assessment  analysis  conducted  on  the  Lerfald  data 
base  is  that  when  dealing  with  nearly  complete,  uncomplicated  and  easy-to-interpret  ionograms,  the 
accuracy  of  the  mass-produced  profiles  is  essentially  the  same  as  our  hand  method  when  processing 
assumptions  were  identical.  However,  the  validity  of  some  of  the  assumptions  is  questionable.  For 
example,  it  has  been  amply  proved  by  rocket  and  incoherent  scatter  measurements  that  considerable 
negative  electron  density  gradients  or  valleys  can  exist,  especially  between  the  E  and  F  layers  during 
deep  night  periods.  When  an  E-layer  model  and  a  monotonic  analysis  are  used  for  these  periods,  the 
resulting  "true"  height  calculation  is  too  low,  and  the  electron  density  gradient  is  too  high.  The 
magnitude  of  the  error  depends  upon  the  depth  and  width  of  the  valley  at  the  time  of  the  ionogram. 
The  resulting  Lerfald  profiles  will  be  the  lowest  practical  approximation  of  the  "true"  heights.  This 
is  not  to  say  they  are  not  useful,  especially  if  the  upper  limit  is  also  known.  The  upper  limit  can 
be  calculated  by  assuming  zero  ionization  below  the  lowest  F  region  echoes  (h'F)  seen  on  the  ionogram. 
In  fact,  the  system  presently  used  by  EDS  calculates  an  upper  and  a  lower  limit  profile  when  there  is 
no  clear  choice  of  a  monotonic  or  nonmonotonic  analysis.  Both  profiles  are   sent  to  the  user,  so  he 
will  know  the  range  of  the  possible  errors. 

A  modified  Lerfald  procedure  has  been  adopted  by  EDS  for  routine  production  of  electron  density 
profiles,  for  periods  when  the  ionosphere  is  considered  monotonic.  When  the  ionosphere  is  considered 
nonmonotonic,  the  alternate  method  mentioned  above  should  be  used  to  calculate  the  profile.  At  times, 
generally  during  sunrise  and  sunset  periods,  when  it  is  questionable  which  process  will  yield  more 
accurate  results,  ionograms  should  be  processed  by  both  methods,  and  both  results  should  be  considered 
by  the  user.  This  would  give  him  a  calculation  for  the  highest  and  lowest  reasonable  profile.  It  has 
been  suggested  that  all  profiles  should  be  calculated  this  way.  The  value  of  the  Lerfald  data  base 
would  be  increased  many  times  over  if  each  ionogram  had  been  processed  for  both  the  upper  and  lower 
limit  profiles. 
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APPENDIX 
Noon  and  Midnight  Excerpts  from  Data  Base 


This  Appendix  displays  all  of  the  noon  and  midnight  (local  time)  profiles  resulting  from  the 
Lerfald  analyses  described  in  this  Report.  The  format  of  the  display  is  as  shown  in  Figure  1.  For 
easy  use  of  this  Appendix,  an  index  is  provided  below. 


Station  Dates  -  1964 

BC  Boulder       Mar  17,  18,  19 

Jun  16 

Jun  17,  18,  19  (noon  only) 

Sep  20 

Sep  21 

Sep  22,  23  (midnight  only) 

Sep  24 

Dec  10 

Dec  11,  12  (midnight  only) 

Dec  13 

Dec  14 

Dec  15  (midnight  only) 

Dec  16  (noon  only] 

Dec  17 


Page 

Station 

23 

AD 

Adak 

24 

25 

26 

27 

Dates  -  1964 


Page 


OK 

Okinawa 

Mar 

17 

(noon  only) 

28 

Mar 

18, 

19 

PA 

Point  Ar 

Mar 

20 

(midnight  only) 

Jun 

16 

(noon  only) 

29 

Sep 

20 

Sep 

21, 

22  (noon  only) 

Sep 

23 

(midnight  only) 

30 

Sep 

24 

Dec 

10 

(noon  only) 

Dec 

11 

(midnight  only) 

Dec 

13, 

14,  15,  16 

31 

Dec 

17 

32 

HU 

Huancayo 

Dec 

18 

(midni  ght  only) 

GB 

Grand  Bahama 

Mar 
Jun 
Jun 
Sep 
Sep 
Dec 
Dec 

17, 

16 

17 

19, 

22 

10, 

13 

18,  19 

20,  21 
11,  12 

33 
34 
35 
36 

Dec 

14 

(noon  only) 

MA 

Maui 

Dec 

15, 

16 

Dec 

17 

37 

WS 

White  Sands 

Mar 
Mar 
Jun 
Jun 

17 

18 
17 

18 

(midnight  only) 
(noon  only) 

38 

Sep 
Sep 
Sep 

19, 

21 

22 

20 

(noon  only) 
(midnight  only) 

39 

Sep 

23, 

24 

40 

Dec 

10 

(midnight  only) 

Dec 

11 

Dec 

12, 

13,  14 

41 

Dec 

15 

(noon  only) 

Dec 

16, 

17 

42 

Mar  16 

(noon  only) 

43 

Mar  17, 

18 

Jun  15 

(noon  only) 

44 

Jun  16, 

17 

Jun  18 

(noon  only) 

Sep  18 

(noon  only) 

45 

Sep  19 

Sep  20 

(noon  only) 

Sep  21 

Sep  22, 

23,  24 

46 

Dec  9 

(noon  only) 

47 

Dec  10, 

11,  12 

Dec  13, 

14,  15  (noon  only) 

48 

Dec  16 

Dec  17  49 

Mar  17,  19  (midnight  only)   50 

Jun  16  (noon  only) 

Jun  17  (midnight  only) 

Jun  18  (midnight  only)     51 

Sep  19,  20,  21 

(midnight  only) 
Sep  22,  23  (midnight  only)  52 
Dec  12,  13 
Dec  16,  17  (midnight  only)  53 

Jun  16,  17  (midnight  only)  54 

Jun  18 

Sep  23  (midnight  only) 

Sep  24  (midnight  only)     55 

Dec  10,  11,  12  (noon  only) 

Dec  13  (midnight  only)     56 

Dec  14,  15 

Dec  16  (noon  only) 

Dec  17  (noon  only)        57 

Mar  17,  18,  19  58 

Jun  16 

Jun  17  (noon  only)        59 

Jun  18  (midnight  only) 

Sep  19  (midnight  only) 

Sep  20  (noon  only) 

Sep  22  60 

Sep  23  (noon  only) 

Dec  10,  11 

Dec  13  (midnight  only)     61 

Dec  14 

Dec  15,  16  (noon  only) 

Dec  17  62 
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Station 


Dates  -  1964 


Pac 


Station 


CP  Concepcion  Mar  17  (noon  only) 

Mar  18 

Mar  19  (noon  only) 

Jun  16 

Jun  17 

Sep  19  (midnight  only) 

Sep  20,  22  (noon  only) 

Sep  23  (noon  only) 

Sep  24 

Dec  10  (midnight  only) 

Dec  11 

Dec  12  (midnight  only) 

Dec  15,  16,  17  (noon  only) 

MN  Manila     Mar  17  (noon  only) 
Mar  18 

Mar  19  (noon  only) 
Jun  17 

Jun  18  (noon  only) 
Sep  20 

Sep  21,  22  (noon  only) 
Sep  23 

Sep  24  (noon  only) 
Dec  10  (noon  only) 
Dec  11 

Dec  12,  13,  14,  15 
Dec  16 
Dec  18  (noon  only) 


64 

65 

66 

67 


Dates  - 

1964 

Mar 

17, 

18 

Mar 

19 

(noon 

only) 

Jun 

16 

Jun 

17, 

18 

Sep 

19 

(midn 

ght  only 

Sep 

22 

(noon 

only) 

Sep 

23 

Sep 

24 

(noon 

only) 

Dec 

10 

(noon 

only) 

Dec 

11 

Dec 

12 

Dec 

13, 

14  ( 

ioon  only 

Dec 

15 

Dec 

16 

Dec 

17 

(noon 

only) 

Page 
72 

73 

74 

75 


69 


70 

71 


JA  Jamaica  Jun  16,  17,  18  (midnight  only) 

Sep  19 

Sep  20,  21,  22,  23  (midnight 

only) 

Sep  24  (noon  only) 

Dec  10  (midnight  only) 

Dec  11  (noon  only) 

Dec  12 

Dec  14  (noon  only) 

Dec  15  (midnight  only) 


77 
78 
79 


22 


EXCERPTS   FROM  DATA  BASE 
NOON  MIDNIGHT 


400 


-  soo 


X   200 

0 
3 

*■   100 


1/ 


if    id1 

ELECTRONS  PC*  CC 


10' 


\t         to' 

HAOC-  I 

(IC  840317  1201) 


400 


'« 

X 

® 

3 


10*      1  o'      10' 

ELECTRONS  PER  CC 


10 


10*     10' 

SRM>E»  1 

(IC  C40319  1200) 


\t  \t  18* 

ELECTRONS  PER  CC 


10 


II* 
6JU0t>  1 

(IC  C40SIS  1200) 


10*      II* 
ELECTRONS  PER  CC 


«UDE»  9 

(IC  S40C1C  1200) 


400 


400 


300 


200 


ELECTRONS  PER  CC 


\i  10*       10* 

ELECTRONS  PER  CC 


tO       10 
ELECTRONS  PER  CC 


10       10 
ELECTRONS  PER  CC 


10*      10' 
«AOE«  1 

IIC  (40317    0! 


10 


6RA0C-  I 

(IC  C403II    0) 


10       10      10 
CRAPE*  1 
IIC  C40319    0) 


10     to 

SSAOE'  1 

IIC  C40C1C    0) 
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EXCERPTS  FROM  DATA  BASE 


\i 


NOON 


io*     \t  to'    to' 

MADE-  1 
ELECTRONS  PER  CC  (DC  S40C17  Ittll 
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300 

200 

/ 

I 

5 

too 

J 

1 

i 

if 

10* 

10* 

10* 

10* 

10       10 
ELECTRONS  RER  CC 


10 


MADE-  t 

ELECTRONS  RER  CC  (»C  (40C1I   1100) 


10'  10 

WADE-   1 
IIC  (40919  1200) 


400 

1 

soo 

200 

J 

% 

w 
1- 

100 

•— 

J 

0* 

\$ 

10* 

10' 

10* 

10' 

ELECTRONS  RER  CC 


MADE"   I 

IIC  840920    1200) 


MIDNIGHT 


\i  10*  10* 

ELECTRONS  RER  CC 


10*  10*  10 

6RA0E«  2 
IIC  (40920         01 
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400 


\i 


3 


I 
3 


EXCERPTS   FROM  DATA  BASE 
NOON  MIDNIGHT 

400 
100 


200 


lo"     10" 
electrons  per  cc 


10 


10       10" 
ELECTRONS  PER  CC 


10     to' 

6RA0E*  t 

IIC  (40921  1200) 


10       10"      10' 
WADE-  1 
<IC  (40924  1200) 


10'      to 
ELECTRONS  PER  CC 


10      10 
IIC  (40921    0) 


400 

301 

2(1 

100 

1 

i 

\i 

ti* 

10* 

\i 

10' 

WW  1 

ELECTRONS  PER  CC  <IC  141111    0) 


400 


10*      101      10* 
ELECTRONS  PER  CC 


io'     to*    10* 

MADE"  I 

IIC  (40925    0) 


400 


300 


10* 


\i 


ELECTRONS  PER  CC 


10*      10 
(ROE*  2 
IIC  (41924    0) 
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3 


«> 

I 


I 
a> 


« 

I 


EXCERPTS   FROM  DATA  BASE 
NOON  MIDNIGHT 


l»"  il* 

ELECTRONS  »ER  cc 


10 


10*  10*  10* 

ELECTRONS  PER  CC 


10 


ti*     10' 

HAOC"  t 

<IC  C412I0  ttlOl 


10      10 
«UDE»  1 
t|C  C4I21S  1200) 


\t  \t  \< 

ELECTRONS  n%  cc 


to*         to' 

tIC  C41I10         I) 


411 

311 

2*1 

til 

I 

i 

\i 

\< 

\i 

if 

ti» 

ELECTRONS  n*  cc 


MAM*  I 

tIC  Mlltt         I) 


411 

J00 

211 

til 

\i  ti*  ti*  i«T  ti*         to' 

MADE*  t 
ELECTRONS  »ER  CC  tIC  C4I2I2    I) 


411 

Sll 

211 

111 

1 

t 

\t 

If* 

to* 

to* 

ll' 

ELECTRONS  »ER  CC 


MAOC'  I 

tIC  C4I21S    0) 
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EXCERPTS   FROM  DATA  BASE 
NOON  MIDNIGHT 


10'     to' 

10' 

10"      10' 

«AD€«  1 

\i 

it               10* 

10* 

to1 

6RADE"  1 

ELECTRONS  PER  CC 

IIC  C412U  1200) 

ELECTRONS  RE*  CC 

(IC  (41214 

if  10'  10" 

ELECTRONS  RER  CC 


10 


10  10' 

MADE*   1 

IIC  (41216   1200) 


10  10' 

ELECTRONS  RER  CC 


10 


10  10' 

UAOE>   1 

IIC  (41217   1200) 


400 


300 


10*      \t 


ELECTRONS  RER  CC 


10 


1 1*     to' 
6RAOC*  I 

(IC  (41211    0) 


ELECTRONS  RfR  CC 


II*  II- 

WAN>  t 

IIC  (41217         0) 


27 


EXCERPTS  FROM  DATA  BASE 
NOON  MIDNIGHT 


1/  if  l/  if 

ClCCTRMS  PfM  CC 


1/  If 

•  I 
(OK  MIII7  Kill 


\i  \t  if  if  if  if 

MAM"  1 
ELECTRONS  Hk  CC  COR  MIHI  HID 


if 


10  10" 

ELECTRONS  "ER  CC 


II*  II' 

MADE"  I 
(OK  C405IJ  1211) 


I 


411 
311 
211 

III 


i    i  i  i  ni«       i    i  r r> 


1/  if  if  if  »f  »•' 

MAM"  I 

CUCTKOMI  »t»  CC  IOK  (41119         I) 


411 

311 

211 

III 

1 

f 

II* 

1/ 

II* 

II* 

II* 

ELECTRON!  »ER  CC 


MAN*  I 

(OK  C4I3II         I) 


411 

311 

211 

III 

t 

f 

if 

II* 

if 

ii* 

if 

MAM"  I 

ELECTRONS  PER  CC  IOK  C4I320         0> 
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EXCERPTS  FROM  DATA  BASE 


NOON 

400 

*.    100 

o 
2   too 

/ 

> 

J 

^     100 

S 

MIDNIGHT 


io*  to* 

ELECTRONS  PER  cc 


10'     10*     to* 
MADE"  I 
(OK  I40IIC  1100) 


10       10 
ELECTRONS  PER  CC 


to 


10       10 
ELECTRONS  PER  CC 


10      10 
6RA0E«  1 
(OK  (40920  1200) 


10       10      10 
«UDE»  1 

(OK  $49921  1200) 


ELECTRONS  PER  CC 


10*    to' 
MADE"  1 
(OK  (40122  1200) 


ELECTRONS  PER  CC 


to    to 

CRAOE*  I 

OK  (40120    0) 


29 


EXCERPTS   FROM  DATA  BASE 


NOON 


MIDNIGHT 


I 


I 

3 


10'      10* 
ELECTRONS  »ER  CC 


\t  to* 

MAN'  I 

(OK  (40114  1I0II 


400 

300 

5  200 

J 

> 

5 

p  too 

J 

i 

to' 

to* 

to' 

to1 

10' 

ELECTRONS  »ER  CC 


WOE*  I 

(OK  C4I2I0   ttOO) 


to*  \i  to*  io*  to*         \i 

MAOS*  t 

ELECTRONS  »ER  CC  IOK  MOMS         0) 


ELECTRONS  »ER  CC 


(RASE*   1 

(OK  C40S44         0) 


400 

300 

200 

too 

1 

i 

\$ 

\( 

to' 

to' 

to' 

ELECTRONS  *ER  CC 


MADE1   t 

(OK  (4121!  0) 
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EXCERPTS   FROM  DATA  BASE 


NOON 


400 


«.      300 


X      200 


100 


1/ 


a     a    i  i  him         r     i    -  *•***■    m 

io'  id1 


ELECTRONS  PER  CC 


10" 


to*     11* 
MADE*  I 

(OK  141213  11001 


400 

300 

^ 

9 
€> 
X 

200 

100 

■ 

1 

1 

\i 

\i 

10* 

\t 

io* 

ELECTRONS  PER  CC 


SRAOS-  1 

(OK  (41214  1200) 


400 

300 

9 

) 

X 

200 

/ 

r 

3 

100 

y 

1 

i 

io' 

10* 

\t 

\i 

\i 

ELECTRONS  PER  CC 


HAOE*  1 

(OK  C4I21I  HID 


400 


X 
a> 

3 


io1  to'  10* 

ELECTRONS  PER  CC 


10 


\i         io' 

6RADE-  2 

(OK  S4I2IC  1209) 


MIDNIGHT 


490 
3t» 

200 

100 


\i 


10 


ELECTR0W3  PER  CC 


If     10 
S«40E»  t 
(OK  841213    0) 


10'      10* 
ELECTRONS  PER  CC 


o1 

«AOf«  1 

(OK  841214         0) 


400 

sot 
1(1 

i  eof 


is? 


r     i  unm        n    nimi        *    *-  •  ■'~jl 


i  §r 


is- 


is 


ELECTRONS  PER  CC 


\i        io' 

(RACE"  1 

(OK  S412IB    01 


10*      10*      10* 
ELECTRONS  PER  CC 


10 


10*  10 

WAOE"  1 

(OK  C4I21S         0) 
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EXCERPTS   FROM  DATA  BASE 


NOON 


MIDNIGHT 


\t 


ELECTRON!  »ER  CC 


to'  \t  II' 

MADE"  I 

(OR  I4IIIT   till! 


5 

I 


I 


X 

3 


\i  \f  \i  \t  \t        u 

•MAN*  I 

ELECTRON!  MR  CC  I0R  Mill  7         I) 


in       i    ilium       ii    ill  inni        i 

\i  \i  \t  \i  \t        to' 

ELECTRONS  MR  CC  (OR  Ml  til         t> 
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1/ 


EXCERPTS   FROM  DATA  BASE 
NOON  MIDNIGHT 


400 

*- 

300 

5 

X 

200 

100 

1 

t 

.0' 

10* 

to' 

10* 

10' 

ELECTRONS  PER  CC 


6RA0E»  1 

((•  (4031?  12001 


400 

300 

♦- 

S 

I 

200 

® 

3 
w 

100 

1 

t 

to1 

10' 

10* 

1/ 

10' 

ELECTRONS  PER  CC 


MAOC-  t 

((•  84J31I  12101 


10'     II*     to' 
(*ADC»  I 

electrons  per  cc  tsi  uostt  1200 


10"      Iff 
ELECTRONS  PER  CC 


.1* 


II*     !•' 
«RAK>  I 
((»  (40(l(  till) 


400 

300 

200 

100 

1 

i 

\i 

10' 

10* 

10* 

10* 

ELECTRONS  PER  CC 


400 


\i  10'       10' 

ELECTRONS  PER  CC 


ELECTRONS  PER  CC 


400 


300 


211 


100 


(RU>E»  1 

((•  (40317    0) 


(RACE*  1 

ft*  (40311    0) 


400 

300 

200 

100 

1 

i 

to' 

10* 

.0* 

10* 

!•' 

6RA0C*  I 

(«•  (4*319    I) 


1/ 


II*       10* 
ELECTRONS  PER  CC 


it 


\t  II* 

(ROE*  I 

<(*  (40(1 (    I) 
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EXCERPTS   FROM  DATA  BASE 


NOON 


MIDNIGHT 


\i  \t  io* 

ELECTRONS  m   CC 


X 


II 


1/      it 

6«ADf>  I 

I6fi  64161?  18*1) 


10       10 
ELECTRONS  PER  CC 


10 


II*     10* 
6RA0S*  I 
(61  64(119  120D 


40C 

300 

•*- 
£ 

,2" 
'5 

a: 

200 

^y 

; 

© 

3 

1  00 

^-^ 

1 

0* 

to' 

10* 
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«/ 

10* 
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GRADE*  1 
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1 
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/ 

) 
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*• 

^ 

1 

i 

\i 

19* 

II* 

it 
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SRAOC*  I 

(6§  641921  till) 


411 
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CRAM*  I 
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it 

it 
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6RA0C*  1 

(61  641921         I) 
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EXCERPTS  FROM  DATA  BASE 
NOON  MIDNIGHT 
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X 
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10*  10*  \t  io' 
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NOON 

5 

X 

400 

soo 

200 

j 

; 
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to' 
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ELECTRONS  PER  CC 
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10 


to     to" 
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to 


«ADC»  t 
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\t        to- 

SRAM-  I 
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«AO€«   t 

<(•  841211   1201) 
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soo 
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/ 

; 
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w. 

t- 

y 

1 

<f 
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to' 

to1 

to1 
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ELECTRONS  RER  CC 


«AOC»  1 
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if 
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l/ 
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\t 

\i 
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1/ 
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NOON  MIDNIGHT 


10'  10 
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6RA0C"  I 
(OS  C40930  11001 


400 

SOO 

Height 

o 

y 

True 

O 
IB 

y 

10'      10* 
ELECTRONS  PER  CC 


10 


10      10 
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MADE*  1 

INS  C40920    0) 


40C 

300 

200 

100 

1 

i 

\i 

10* 

.0' 

10* 

to' 

CRAOC"  I 

IVS  C409S2    0) 


39 


EXCERPTS   FROM  DATA  BASE 


NOON 


MIDNIGHT 


if  if  if  l/  if  if 

MAN*  I 
ELECTRONS  RIR  CC  l»t  MINI  ltll> 


400 

300 

1 

I 

200 

I 

1 

100 

y 

1 

f 

if 

ii* 

\i 

if 

11' 

ELECTRONS  RER  CC 


X 
■3 


It       10 
ELECTRONS  PER  CC 


10 


«AOf  •   I 

(MS  Ml«4   till) 


10  10" 
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NOAA  SCIENTIFIC  AND  TECHNICAL  PUBLICATIONS 

NO  A  A,  the  National  Oceanic  and  Atmospheric  Administration,  was  established  as  part  of  the  Depart- 
ment of  Commerce  on  October  3,  1970.  The  mission  responsibilities  of  NOAA  are  to  monitor  and  predict  the 
state  of  the  solid  Earth,  the  oceans  and  their  living  resources,  the  atmosphere,  and  the  space  environment  of 
the  Earth,  and  .to  assess  the  socioeconomic  impact  of  natural  and  technological  changes  in  the  environment. 

The  six  Major  Line  Components  of  NOAA  regularly  produce  various  types  of  scientific  and  technical 
information  in  the  following  kinds  of  publications: 


PROFESSIONAL  PAPERS  —  Important  definitive 
research  results,  major  techniques,  and  special  in- 
vestigations. 

TECHNICAL  REPORTS— Journal  quality  with  ex- 
tensive details,  mathematical  developments,  or  data 
listings. 

TECHNICAL  MEMORANDUMS  —  Reports  of 
preliminary,  partial,  or  negative  research  or  tech- 
nology results,  interim  instructions,  and  the  like. 

CONTRACT  AND  GRANT  REPORTS— Reports 
prepared  by  contractors  or  grantees  under  NOAA 
sponsorship. 


TECHNICAL  SERVICE  PUBLICATIONS— These 

are  publications  containing  data,  observations,  in- 
structions, etc.  A  partial  listing:  Data  serials;  Pre- 
diction and  outlook  periodicals;  Technical  manuals, 
training  papers,  planning  reports,  and  information 
serials;  and  Miscellaneous  technical  publications. 


ATLAS — Analysed  data  generally  presented  in  the 
form  of  maps  showing  distribution  of  rainfall,  chem- 
ical and  physical  conditions  of  oceans  and  atmos- 
phere, distribution  of  fishes  and  marine  mammals, 
ionospheric  conditions,  etc. 
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Information  on  availability  of  NOAA  publications  can  be  obtained  from: 

ENVIRONMENTAL   SCIENCE   INFORMATION   CENTER 

ENVIRONMENTAL   DATA  SERVICE 

NATIONAL  OCEANIC  AND  ATMOSPHERIC  ADMINISTRATION 

U.S.   DEPARTMENT  OF   COMMERCE 

3300  Whitehaven  Street,  N.W. 
Washington,  D.C.  20235 


